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In  this  work,  we  present  the  Poroelastic  Relaxation  Indentation  (PRI)  testing  approach  for  quantifying 
the  mechanical  and  transport  properties  of  thin  layers  of  poly(ethylene  glycol)  hydrogels  with 
thicknesses  on  the  order  of  200  pm.  Specifically,  PRI  characterizes  poroelastic  relaxation  in  hydrogels 
by  indenting  the  material  at  fixed  depth  and  measuring  the  contact  area-dependent  load  relaxation 
process  as  a  function  of  time.  With  the  aid  of  a  linear  poroelastic  theory  developed  for  thin  or 
geometrically  confined  swollen  polymer  networks,  we  demonstrate  that  PRI  can  quantify  the  water 
diffusion  coefficient,  shear  modulus  and  average  pore  size  of  the  hydrogel  layer.  This  approach 
provides  a  simple  methodology  to  quantify  the  material  properties  of  thin  swollen  polymer  networks 
relevant  to  transport  phenomena. 

Introduction 

Thin  layers  of  swollen  polymer  networks  are  ubiquitous  in 
various  technological  applications  relevant  to  sustainability, 
energy,  and  bioengineering.  Some  examples  include  membranes 
for  water  purification,1,2  proton  exchange  membranes  used  in 
fuel  cells,3,4  scaffolds  for  tissue  engineering,5  and  polymer 
particles  for  drug  delivery.6,7  In  these  examples,  their  perfor¬ 
mance  is  closely  linked  to  the  ability  to  regulate  transport  of 
small  molecules  such  as  solvent.  Hence,  understanding  the 
critical  time-  and  length-scales  that  regulate  transport  of  these 
molecules  within  the  polymer  network  provides  vital  informa¬ 
tion  to  enable  development  of  membrane  materials  with 
improved  performance. 

Various  techniques  including  impedence-,8,9  optical-,10,11 
scattering-,12  and  spectroscopic-based13-16  measurement 
approaches  have  been  used  to  quantify  the  through-thickness 
solvent  diffusion  within  swollen  polymer  network  layers. 

Since  the  characteristic  solvent  diffusion  time  scales  with  the 
film  thickness  squared,  this  implies  that  through-thickness 
diffusion  measurements  become  more  challenging  for  thinner 
films. 

Recently,  indentation  testing  have  been  demonstrated  as  an 
interesting  approach  to  quantify  solvent  transport  in  polymer 
gels.17-23  An  illustrative  example  is  presented  in  Fig.  1.  The 
approach  involves  indenting  a  hydrogel  with  a  spherical  probe  to 
a  fixed  depth  (5)  (Fig.  la).  In  response  to  this  indentation, 
a  corresponding  load  develops  that  is  defined  by  the  hydrogel 
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stiffness.  As  testing  time  progresses,  the  hydrogel  will  undergo 
a  stress-relaxation  process  (Fig.  lb).  Generally,  two  primary 
relaxation  processes  can  occur  that  are  related  to  the  dynamics 
on  the  polymer  chain  level.  The  first  process  is  viscoelastic 
relaxation,  and  is  associated  with  conformational  changes  of  the 
polymer  chain  due  to  the  applied  deformation.  The  second 


solvent  reservoir 


Fig.  1  Spherical  indentation  testing  of  relaxation  phenomena  in 
hydrogels,  (a)  A  spherical  probe  is  brought  into  contact  with  the  hydrogel 
surface  in  a  water  environment.  The  probe-hydrogel  and  hydrogel- 
substrate  interfaces  are  considered  frictionless.  The  spherical  volume 
directly  beneath  the  probe  represents  the  hydrogel  region  where  solvent 
migration  occurs,  (b)  Compressing  the  gel  at  a  fixed  displacement  (5) 
results  in  a  relaxation  in  load  (P)  that  is  characterized  by  a  relaxation 
time  (r). 
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process  is  poroelastic  relaxation.  This  relaxation  process  is 
associated  with  the  transport  of  water  molecules  out  of  the  region 
of  hydrogel  material  under  deformation  as  defined  by  the  contact 
area  Tt a1.  Both  processes  are  concurrent,  with  each  leading  to 
a  macroscopic  relaxation  of  the  load  response.  However,  the 
characteristic  relaxation  time-scale  (t)  of  each  process  is  unique 
since  the  two  relaxation  times  depend  on  different  sets  of  mate¬ 
rials  properties. 

The  viscoelastic  relaxation  time-scale  {iR)  is  determined  by  the 
ratio  of  the  shear  viscosity  (77)  and  shear  modulus  (G)  of  the 
hydrogel,  iR  ~  77/G.24  Both  of  these  material  properties  are 
independent  of  the  contact  radius  established  by  a  typical  milli¬ 
metre-sized  indenter  ( a  ~  100  pm).  This  implies  that  iR  is  also 
independent  of  the  contact  size. 

On  the  other  hand,  the  poroelastic  relaxation  time-scale  (iP)  is 
proportional  to  the  size  of  the  contact.  This  contact-size  depen¬ 
dence  is  attributed  to  the  time-scale  required  for  the  solvent 
within  the  deformed  volume  of  hydrogel  material  to  migrate  over 
a  distance  defined  by  the  contact  radius.  Specifically,  zP  is  related 
to  the  contact  radius  and  the  diffusion  coefficient  (D)  of  the 
solvent  within  the  hydrogel  by:18,19 


Therefore,  the  time-scales  of  the  two  relaxation  processes  are 
distinguishable  based  on  the  contact  size.  As  an  example,  for 
a  typical  swollen  polymer  network  with  parameters  of  D  =  1CF10 
m2/s  and  a  =  100  pm,  the  poroelastic  relaxation  time  (zP  ~  100  s) 
is  significantly  longer  than  its  viscoelastic  relaxation  time  (zR  ~ 
0.1  s).  More  importantly,  eqn  (1)  suggests  that  we  can  adjust  the 
poroelastic  relaxation  time  simply  by  controlling  the  contact 
area. 

This  indentation  testing  approach  is  attractive  for  measuring 
the  material  properties  of  thin  swollen  polymer  networks  for 
several  reasons.  First,  the  test  measures  solvent  transport  over 
the  lateral  dimensions  of  the  polymer  film  as  opposed  to  the 
through-thickness  dimension.  For  polymer  films  with  isotropic 
properties,  this  measurement  approach  is  not  limited  by  the 
film  thickness  of  the  polymer  layer  as  compared  to  the  other 
techniques.  Second,  the  approach  is  quite  simple  and  straight¬ 
forward  to  implement  since  it  only  involves  indenting  the 
material  at  fixed  depth  with  a  macroscopic  spherical  probe. 
Therefore,  the  test  facilitates  direct  measurements  of  the 
materials  properties  of  the  polymer  film  that  is  supported  by 
a  substrate.  Third,  several  material  properties  including  diffu¬ 
sion  coefficient,  shear  modulus,  and  average  pore  dimension 
can  be  extracted  directly.  We  envision  adapting  this  technique 
for  membranes  because  it  can  directly  measure  the  poroelastic 
relaxations  of  a  membrane  film  and  correlate  the  response  to 
material  properties  such  as  diffusivity  and  average  pore 
dimension. 

Although  the  approach  has  been  demonstrated  for  thick 
hydrogel  materials,  its  relevance  for  thin  hydrogel  layer  and 
other  hydrated  polymer  membranes  have  not  been  explored. 
Therefore,  the  objective  of  this  work  is  to  demonstrate 
Poroelastic  Relaxation  Indentation  (PRI)  in  quantifying  the 
transport  properties  for  thin  or  geometrically  confined  swollen 
polymer  network  based  on  polyethylene  glycol)  (PEG) 
hydrogel  layers. 


Materials  and  methods! 

Mechanics  of  indenting  a  thin  poroelastic  layer 


Indentation  of  a  highly  swollen  polymer  network,  such  as 
a  hydrogel,  is  unique  in  that  the  mechanical  deformation  is 
coupled  to  the  mass  transport  of  solvent  through  the  network, 
which  is  a  phenomena  known  as  poroelasticity.25  This  defor¬ 
mation  causes  a  chemical  potential  gradient  to  develop  within  the 
hydrogel.  Over  time,  the  hydrogel  establishes  a  new  chemical 
equilibrium  via  migration  of  water  away  from  the  region  under 
deformation.  Due  to  this  diffusion  process,  the  hydrogel 
undergoes  a  load  relaxation  that  reflects  the  time-scale  to 
establish  this  new  equilibrium,  which  corresponds  to  the  poroe¬ 
lastic  relaxation  time. 

Immediately  following  indentation,  the  water  molecules  do  not 
have  sufficient  time  to  migrate  away  from  the  deformed  region  of 
the  hydrogel.  Therefore,  the  hydrogel  is  considered  to  be 
incompressible  (v  =  0.5),  and  the  indentation  response  is  equiv¬ 
alent  to  deformation  of  an  elastic  layer.  For  an  elastic  film 
deformed  by  a  spherical  probe  to  an  indentation  depth  5,  the 
resultant  load  is  defined  as:26 

P=P„-fP{yR5/h )  =  ■  (2) 

where  G  and  v  correspond  to  the  shear  modulus  and  Poisson’s 
ratio  of  the  elastic  material,  respectively.  The  parameter  PH  is  the 
Hertz  load  predicted  by  Hertzian  contact  mechanics  of  an  elastic 
layer  considered  to  be  semi-infinitely  thick  in  relation  to  the 
Hertz  contact  radius  defined  as, 

aH  =  VRS  .  (3) 


Eqn  (3)  does  not  provide  an  accurate  description  of  the  contact 
radius  when  it  becomes  commensurate  or  greater  than  the 
hydrogel  layer  thickness.  Here,  the  contact  radius  is  related  to  the 
indentation  depth  by: 


:  =  VRd-fa(VR8/h 


(4) 


The  function  fayyR5/hj  accounts  for  deviations  from  Hertz 
mechanics  where  the  substrate  stiffness  influences  the  geometric 
relationship  between  contact  radius  and  indentation  depth.  This 
function  was  solved  numerically  based  on  the  work  of  Yu  and  co¬ 
workers.26  Details  of  this  calculation  are  discussed  in  a  separate 
contribution.27  The  approximate  analytical  form  is: 


\M(y/M/h)2+0.5l(VR8/h)  +0.5 
(VRt/h)2+0A9(VM/h)  +0.5 


(5) 


Similarly,  the  function  fpyjR5/hj  accounts  for  the  deviation 
from  Hertz  load  when  the  stress  field  is  affected  by  substrate 
stiffness.  We  solve  this  function  numerically  and  details  of  the 
numerical  calculation  are  discussed  separately.27  The  approxi¬ 
mate  analytical  form  is: 


f  The  equipment  and  instruments  or  materials  are  identified  in  the  paper 
in  order  to  adequately  specify  the  experimental  details.  Such 
identification  does  not  imply  recommendation  by  NIST,  nor  does  it 
imply  the  materials  are  necessarily  the  best  available  for  the  purpose. 
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/  ^  \  _  2.36(VR6/h)2+  0.82(^R5/h)  +  0.46 
fp[VR6/  )  -  (VR5/h)  +0.46 


(6) 


Hence,  the  instantaneous  load  response  (P,)  of  the  hydrogel 
immediately  following  indentation  is, 

Pi  =i^Rl/25i/2G-fp(VRS/l^  .  (7) 

As  time  progresses,  the  solvent  molecules  migrate  from 
the  deformed  hydrogel  region  until  a  new  chemical  equilib¬ 
rium  is  reached.  Due  to  this  transport  process,  the  hydrogel 
relaxes  to  an  equilibrium  or  long-time  load  ( Pf )  defined  by 
eqn  (2). 

To  model  the  load  relaxation  process,  we  perform  finite 
element  calculations  using  ABAQUS  software.  In  the  calcula¬ 
tion,  the  probe  indents  the  hydrogel  layer  rapidly,  with 
a  loading  time  10~5  times  faster  than  the  total  relaxation  time. 
Both  the  probe-hydrogel  and  hydrogel-substrate  interfaces  are 
considered  to  be  frictionless  and  impermeable.  Within  the 
contact  area,  the  water  flux  is  zero  along  the  probe-hydrogel 
interface.  The  hydrogel  deformation  follows  the  same  config¬ 
uration  of  the  probe.  Outside  the  contact  area,  the  boundaries 
are  traction-free  and  in  diffusive  equilibrium  with  the  external 
water  reservoir.  To  simulate  a  rigid  substrate  with  a  frictionless 
and  impermeable  hydrogel-substrate  interface,  the  indentation 
of  the  hydrogel  layer  in  vertical  direction  is  constrained  and 
zero  flux  is  established.  Details  of  the  calculations  are 
discussed  in  a  separate  contribution.27  The  calculations  were 
performed  to  determine  the  poroelastic  load  relaxation 
behavior  for  a  broad  range  of  y^RS/h  ratios  from  0.11  to  2.98 
(Fig.  2a).  The  results  are  described  in  terms  of  a  dimensionless 
load  relaxation  function  g,  versus  dimesionless  relaxation 
time,  DtIRS: 

''f  -p1  =  s(Dt/R5,  s/Rd/liJ  .  (8) 


The  results  of  Fig.  2a  indicate  that  the  shape  of  the  load 
relaxation  curve  depends  significantly  on  the  degree  of 
geometric  confinement  as  defined  by  the  s/RS/h  ratio.  An 
interesting  feature  is  that  the  poroelastic  relaxation  time 
required  to  relax  from  P,  to  Pf  descreases  with  increasing 
y^RS/h  ratio,  suggesting  that  the  poroelastic  relaxation  time  for 
a  geometrically  confined  hydrogel  of  a  given  thickness 
decreases  with  indentation  depth. 

To  develop  an  analytical  solution  that  describes  the 
poroelastic  relaxation  behavior  for  the  Vlid/h  ratios  explored, 
we  extrapolate  an  approximate  g  function  from  each  curve.  We 
find  the  g  function  takes  the  stretched  exponential  form, 
commonly  used  to  describe  polymer  chain  dynamics 
(Fig.  2b):28'29 

g(^Dt/R6,VR6/hj  =exp(~a(Dt/RSY^j  .  (9) 

The  parameters  a  and  (3  (Fig.  2b  inset)  are  related  to  the 
y^RS/h  ratios  and  represented  by  polynomials  defined  as: 

a  =  1.15  +  0.44(v/P5//?)  +  0.89  (^Rd/h)2 

(10) 

-0.42(v/P5/A)3+0.06(VP5//2)4  , 


Dt/R5 


a(DVR8)p 

Fig.  2  (a)  Finite  element  calculations  of  the  poroelastic  load  relaxation 
function  g,  for  hydrogel  layers  at  different  y^RS/h  ratios.  The  inset  is  an 
enlarge  view  of  the  same  result,  (b)  Collapse  of  the  entire  set  of  g  func¬ 
tions  with  a  stretched  exponential  function.  The  inset  plot  summarizes  the 
a  (unfilled  symbols)  and  0  (filled  symbols)  parameters,  which  are  fitted 
with  polynomial  functions  as  represented  by  eqn  (10)  and  eqn  (11), 
respectively. 


0  =  0.56  +  0.25  (\/R6/li)  +  0.28  (VM/h)2 

-  0.31  (VRS/h)3+  0.1  (^R6/h)4  (H) 

-0.0l(yPd//j)5  . 

Preparation  of  PEG  hydrogel  layers 

The  PEG  hydrogel  layers  were  synthesized  by  photo- 
polymerizing  poly(ethylene  glycol)  methyl  ether  methacrylate 
monomer  (molecular  mass  =  475  g/mol,  PEGMA),  poly( ethylene 
glycol)  dimethacrylate  crosslinker  (molecular  mass  =  875  g/mol, 
PEGDMA)  (Sigma- Aldrich,  St.  Louis,  MO)  with  de-ionized 
water,  and  using  2,2-dimethoxyphenyl  acetophenone  (DMPA)  as 
the  photoinitiator  (Irgacure,  BASF  Corp.,  Florham,  NJ).  All 
chemicals  were  used  as  supplied.  Monomer  formulation  con¬ 
taining  90/10%  by  mass  of  PEGMA/PEGDMA  was  prepared. 
DMPA  was  added  (0.5%  by  mass)  and  allowed  to  dissolve 
completely.  This  formulation  was  then  diluted  50%  by  mass  with 
de-ionized  water  before  polymerization. 

To  polymerize  hydrogel  layers  of  defined  film  thicknesses,  the 
monomer  formulation  was  filled  between  two  glass  substrates 
with  spacers  ranging  between  150  pm  to  2000  pm.  Photo¬ 
polymerization  was  conducted  using  a  Novacure  2000  ultraviolet 
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mercury  arc  lamp  (EXFO,  Mississauga,  ON)  at  a  light  intensity 
of  10  mW/cm2  for  1200  s.  The  samples  were  delaminated  from  the 
glass  substrates  and  immediately  immersed  into  de-ionized 
water.  Samples  swelled  to  an  equilibrium  state  and  remained  in 
23  °C  water  for  48  h  prior  to  testing.  Water  was  periodically 
replaced  to  remove  any  low  molecular  mass  species  which  could 
diffuse  out  of  the  hydrogel.  The  thicknesses  of  the  hydrogels  were 
determined  by  measuring  film  cross  sections  via  optical  micros¬ 
copy  prior  to  testing.  The  thicknesses  of  PEG  hydrogel 
layers  investigated  were  2682  pm  ±  7  pm,  1338  pm  ±  5  pm,  and 
267  pm  ±  7  pm. 

Poroelastic  Relaxation  Indentation  (PRI) 

Next,  we  measure  the  poroelastic  relaxations  of  the  hydrogel 
layers  using  PRI.  This  test  monitors  the  load  relaxation,  P,  of  the 
hydrogel  as  a  function  of  time  due  to  an  indentation  depth,  5, 
with  a  spherical  probe.  Since  the  testing  mode  is  a  displacement- 
controlled  type,  the  indentation  depth,  and  therefore  the  contact 
area  remains  fixed  over  the  entire  test. 

The  PRI  measurement  approach  is  illustrated  in  Fig.  la.  First, 
an  equilibrated  hydrogel  layer  is  transferred  into  a  water  reser¬ 
voir  enclosure,  and  is  allowed  to  dwell  for  1  h  prior  to  testing. 
Next,  a  rigid  glass  probe  with  radius  of  curvature  R  (=5  mm)  is 
brought  into  contact  with  the  hydrogel  surface  at  a  constant 
indentation  rate  (=20  pm/s)  to  a  predefined  indentation  depth 
that  is  controlled  by  a  nanopositioner.  During  the  entire  test,  the 
indentation  depth  and  load  are  monitored  by  a  fiber  optic 
displacement  sensor  and  load  cell,  respectively.  Each  hydrogel 
sample  was  indented  at  several  different  depths.  Depending  on 
the  indentation  depth  and  hydrogel  layer  thickness,  PRI  test  can 
span  from  15  min  to  2  h.  The  test  is  complete  when  the  load 
relaxation  has  reached  a  long-time  steady-state  value. 

Representative  PRI  testing  results  for  the  1338  pm  thick  PEG 
hydrogel  are  shown  in  Fig.  3.  Upon  indenting  the  hydrogel  to 
a  depth  5 ,  an  instantaneous  load  P,  is  reached.  As  testing  time 
progresses,  poroelastic  relaxation  is  observed,  which  is  charac¬ 
terized  by  a  time-dependent  relaxation  of  the  instantaneous  load. 
Eventually,  the  load  reaches  a  steady-state  lower  limit,  which  is 
defined  as  the  long-time  load  Pf.  By  indenting  to  greater  depths, 
the  magnitude  of  the  instantaneous  and  long-time  load,  along 
with  the  poroelastic  relaxation  time,  will  change  as  defined  by 
eqn  (2),  eqn  (7)  and  eqn  (9),  respectively. 

The  water  diffusion  coefficient  is  determined  by  constructing 
a  master  curve  defined  by  eqn  (9).  The  approach  is  identical  to 
the  procedure  used  to  collapse  all  the  g  functions  from  the  finite 
element  calculations  into  a  single  master  curve  (Fig.  2b).  First, 
the  results  from  the  PRI  tests  are  converted  to  dimensionless  load 
g  =  ((P  —  Pf)l(Pi  —  Pf))  versus  adjusted  time  (t/R5)  curves.  Then, 
the  a  and  /?  coefficients  for  each  PRI  test  are  determined  based  on 
the  experimentally  defined  VRS/h  ratio.  Finally,  a  best  fit  of  each 
dimensionless  load  versus  dimensionless  relaxation  time  curve  is 
obtained  using  eqn  (9)  with  the  diffusion  coefficient  D  as  the  lone 
fitting  parameter.  As  shown  in  Fig.  3b,  all  the  resultant  curves 
collapse  onto  a  master  curve. 

Mass  swelling  ratio.  Mass  swelling  ratios  of  the  PEG  hydrogel 
layers  were  determined  from  the  ratio  of  the  equilibrium  wet 
mass  to  dry  mass.  The  wet  mass  of  approximately  of  1  g 


a(DVRS)13 


Fig.  3  (a)  Representative  PRI  testing  results  for  the  1338  pm  thick  PEG 
hydrogel.  The  test  applies  a  constant  deformation  by  indenting  the 
hydrogel  with  a  spherical  probe  at  a  fixed  indentation  depth  (5).  This 
indentation  process  results  in  a  relaxation  in  load  ( P )  as  a  function  of  time 
(t).  The  poroelastic  relaxation  process  is  characterized  as  the  time 
required  by  the  hydrogel  to  relax  from  P,  to  Pf  in  order  to  establish  a  new 
chemical  equilibrium,  (b)  Extrapolation  of  the  water  diffusion  coefficient 
by  normalizing  the  relaxation  load  versus  time  curves  into  dimensionless 
load  versus  dimensionless  time  curves.  This  approach  is  identical  to  the 
approach  shown  in  Fig.  2b.  By  fitting  all  the  curves  with  eqn  (9),  the 
diffusion  coefficient  can  be  determined  for  each  test.  The  inset  plot 
summarizes  the  a  and  /3  coefficients  for  each  test,  which  depend  on  the 
degree  of  geometric  confinement  defined  by  VR5/li. 

equilibrated  PEG  hydrogel  was  measured  in  air.  The  samples 
were  then  dried  for  72  h  under  vacuum,  and  the  dry  mass 
measured.  The  calculated  swelling  ratios  were  invariant  to  the 
hydrogel  thickness. 

Results  and  discussion 

We  used  the  poroelastic  model  to  determine  the  materials 
properties  of  the  thin  hydrogel  layers.  Since  PRI  couples  contact 
mechanics  to  fluid  transport,  both  the  mechanical  and  transport 
properties  of  the  hydrogel  can  be  obtained  with  this  one  testing 
method  to  quantify  the  hydrogel  shear  modulus  and  the  water 
diffusion  coefficient  separately. 

First,  we  quantify  the  shear  modulus  of  the  hydrogel  layers.  At 
instantaneous  time  when  water  cannot  migrate  and  the  hydrogel 
is  considered  to  be  incompressible,  the  PRI  test  is  equivalent  to 
the  contact  mechanical  testing  of  an  elastic  solid.  In  contact 
mechanical  testing,  the  shear  modulus  is  determined  in 
a  straightforward  manner  by  extrapolating  the  value  from  the 
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load-indentation  depth  curve,  thus  a  series  of  instantaneous  load 
versus  indentation  depth  data  points  obtained  from  PRI  is 
equivalent  to  a  load  versus  displacement  curve  from  a  contact 
mechanical  test.30  Here,  the  shear  modulus  of  the  PEG  hydrogels 
is  determined  from  the  instantaneous  load  versus  indentation 
depth  results  of  the  PRI  tests. 

Specifically,  we  used  eqn  (7)  to  determine  the  shear  modulus 
for  all  the  hydrogel  layers  based  on  the  set  of  Ph  5  values 
obtained  from  the  PRI  experiments  (Fig.  4a).  The  average  shear 
modulus  is,  G  =  1.0  x  105  N/m2  ±  1.5  x  104  N/m2,  which  is 
within  the  range  for  similar  types  of  PEG  hydrogels.31,32 

The  water  diffusion  coefficients  are  determined  using  the 
procedure  outlined  in  the  Materials  and  Methods  section.  As 
summarized  in  Fig.  4b,  the  average  water  diffusion  coefficient  is, 
D  =  4.4  x  10-10  m2/s  ±  7.6  x  HP11  m2/s,  which  is  in  good 
agreement  with  reported  values  for  a  crosslinked  PEG  diacrylate/ 
PEG  acrylate  hydrogel  with  similar  water  content  (D  =  4.8  x 
10  10  m2/s).33 

The  material  properties  obtained  from  PRI  testing  can  be 
linked  with  the  molecular  structure  of  the  hydrogel  that  governs 
transport  if  we  consider  the  diffusion  of  water  through  the 
hydrogel  network  is  analogous  to  flow  through  a  pore  with  an 
average  size  fp.34 

In  this  case,  the  average  pore  size  is  proportional  to  the 
hydrodynamic  permeability.  In  relation  to  poromechanics  of 
hydrogels,  the  average  pore  size  is  related  to  the  materials 
constants  that  can  be  experimentally  measured  from  PRI  testing 
( D,G,v ).  Specifically,  the  pore  size  is  defined  as:19 


fp  =  2 


(r,D(  l-2y)\ 

\2  G(l-v)J 


1/2 


(12) 


where  77  is  the  viscosity  of  water. 

The  Poisson’s  ratio  is  also  determined  by  the  PRI  experiments. 
Traditionally,  Poisson’s  ratio  describes  the  compressibility  of 
a  material  due  to  mechanical  deformation.  Here,  the  change  in 


Fig.  4  Summary  of  materials  properties  obtained  by  PRI  testing,  (a) 
Shear  modulus  of  the  PEG  hydrogel  as  a  function  of  the  degree  of 
geometric  confinement  for  all  the  hydrogel  layers,  (b)  Diffusion  coeffi¬ 
cient  of  water  as  a  function  of  the  degree  of  geometric  confinement  for  all 
the  hydrogel  layers.  The  dashed  lines  represent  the  average  modulus  and 
diffusion  coefficient  determined  from  the  PRI  tests.  The  error  bars,  which 
are  smaller  than  the  symbol  size,  represent  one  standard  deviation. 


compressibility  is  due  to  the  fact  that  the  water  content  within  the 
deformed  region  is  not  conserved  during  PRI  testing.  Thus,  this 
change  in  Poisson’s  ratio  should  be  considered  as  a  change  in 
solvent  volume  fraction.  Specifically,  this  parameter  is  related  to 
the  ratio  of  the  instantaneous  and  long-time  loads,  PJPf  = 
2(1  —  v)  as  defined  by  eqn  (2)  and  eqn  (5). 17,18  Rearranging  this 
expression  leads  to  the  following: 

Pf-  PJ2  =  v  Pf  (13) 

We  use  eqn  (13)  to  determine  the  average  value  of  the  Poisson’s 
ratio  by  plotting  Pj —  P,l 2  vs.  Pf  (Fig.  5)  of  all  the  hydrogel  layers. 
From  the  slope  of  the  curve,  we  find  v  =  0.43.  Substituting  this 
value,  along  with  the  other  materials  constants  for  the  hydrogel 
and  water  (77  =  10~3  N- s/m2),35  we  determine  E,P  =  1.5  nm. 

From  the  swelling  ratio  measurements  of  the  hydrogel  layers, 
a  volumetric  swelling  ratio,  S  ~  5,  is  obtained.  Based  on 
a  modified  Flory-Rehner  theory  developed  for  isotropic  swelling 
of  hydrogels,  the  average  pore  size  of  a  polymer  network  is 
estimated  to  be:36 


?  =  (S)m(2CnM J Mr)'n  l  (14) 

where  C„  is  the  Flory  characteristic  ratio,  the  ratio  MJMr  is  the 
number  of  C-C  bonds  between  crosslinks,  and  l  is  the  C-C  bond 
distance.  Assuming  the  following  PEG  polymer  chain  constants, 
C„  =  3.8, 37  MJMr  =  725  g  moF744  g  moF1  =  16,38  and  /=  1.54 
A,  we  find  E  ~  2  nm.  We  note  that  the  hydrogel  network  model 
used  in  eqn  (14)  assumes  that  the  polymer  chains  between 
network  junctions  consist  of  linear  homopolymer  chains, 
whereas  the  PEG  hydrogels  used  in  this  work  consist  of  essen¬ 
tially  co-monomers.  Thus,  eqn  (14)  can  only  be  used  as  an  esti¬ 
mation  of  the  average  pore  size.  Nevertheless,  the  pore  size 
determined  by  both  techniques  are  in  close  agreement,  thus 
demonstrating  that  the  PRI  technique  can  be  used  to  estimate  an 
average  pore  dimension  of  the  PEG  hydrogel. 

Conceptually,  it  is  reasonsable  to  expect  that  the  contact 
radius  defines  the  primary  length-scale  that  determines  the 
poroelastic  relaxation  time  irrespective  of  the  hydrogel  layer 
thickness.  This  concept  can  be  demonstrated  by  comparing  the 
PRI  testing  of  a  thick,  or  geometrically  unconfined  (a  <C  h )  and 
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Fig.  5  Determination  of  the  Poisson’s  ratio  upon  complete  relaxation  of 
the  load.  By  plotting  Pf  —  PJ2  vs.  P^of  all  the  hydrogel  layers,  the  slope  of 
the  curve  defines  the  average  value  of  the  Poisson’s  ratio,  v  =  0.43. 
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a  thin,  or  highly  geometrically  confined  ( a  >  h)  hydrogel  layer 
(Fig.  6). 

In  both  cases,  the  indentation  process  deforms  a  volume  of 
hydrogel  material  ( V)  that  is  related  to  the  size  of  the  contact 
radius.  This  volume  is  approximately  a3  and  V ~  a2h,  for  the 
thick  and  thin  hydrogel  layers,  respectively.  In  response,  the 
hydrogel  layers  relieve  the  deformation  via  water  migration  away 
from  the  region  of  indentation.  Specifically,  a  volume  of  water 
must  migrate  over  a  surface  area  (As)  that  is  also  defined  by  the 
shape  of  the  deformed  region.  This  surface  area  is  approximately 
As  ~  a2  and  As  ~  ah,  for  the  thick  and  thin  hydrogel  layers, 
respectively.  Regardless  of  the  hydrogel  layer  thickness,  the 
volume  defines  the  amount  of  water  that  will  migrate  while  the 
surface  area  defines  the  time-scale  of  water  migration.  The 
volume-to-surface  area  ratio  defines  the  lateral  length-scale  for 
solvent  transport  through  the  hydrogel  during  PRI  testing.  Thus, 
this  primary  length-scale,  as  defined  by  a,  determines  the  average 
characteristic  time-scale  of  the  poroelastic  relaxation  process. 

Finally,  Fig.  6  also  qualitiatively  explains  the  finite-element 
predictions  of  Fig.  2  where  we  find  that  the  average  poroelastic 
relaxation  time  decreases  with  increasing  geometric  confinement. 
More  importantly,  the  (}  coefficient,  which  describes  the  breadth 
of  the  distribution  of  poroelastic  relaxation  time,  approaches 
unity  with  geometric  confinement.  This  trend  is  associated  with 
the  relationship  between  pressure  distribution  and  geometric 
confinement.  According  to  Darcy’s  law,  the  pressure  distribution 
is  directly  proportional  to  water  migration.  Since  the  load 
relaxation  process  measured  by  PRI  provides  an  indirect 
measure  of  transport,  any  change  in  the  pressure  distribution, 
such  as  the  extent  of  geometric  confinement,  will  lead  to  changes 
in  the  poroelastic  relaxation  process.  In  the  geometrically 
unconfined  case,  the  pressure  distribution  is  spherical  due  to  the 
shape  of  the  probe  and  absence  of  substrate  stiffness  effect 
(Fig.  6a).  As  a  result,  water  migration  will  be  different  depending 
on  the  location  of  the  water  molecules  within  the  deformed 
hydrogel  region.  At  the  other  extreme  of  highest  geometric 
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Fig.  6  Comparison  of  the  PRI  testing  of  (a)  semi-infinitely  thick 
hydrogel  where  the  contact  radius  is  significantly  less  than  the  layer 
thickness  (a  -C  /;)  and  (b)  thin  film  case  where  the  contact  radius  is  greater 
than  the  layer  thickness  ( a  >  h).  The  volume  directly  beneath  the  probe 
represents  the  hydrogel  region  where  solvent  migration  occurs. 


confinement,  the  pressure  distribution  has  a  more  cylindrical  or 
uniform  profile  due  to  substrate  stiffness  effect  (Fig.  6b).  Here, 
the  uniform  pressure  distribution  implies  that  water  migration 
will  be  independent  of  location  within  the  deformed  hydrogel. 
The  observation  that  the  exponent  (eqn  (9)),  typically  used  to 
describe  the  distribution  of  relaxation  times  in  polymer 
dynamics,  approaches  unity  as  \[R8/h  increases  supports  this 
concept  that  the  breadth  in  the  distribution  of  water  migration 
times  narrows  with  geometric  confinement. 

This  result  has  important  implications  from  a  measurement 
standpoint.  Unlike  other  techniques  for  measuring  diffusion, 
PRI  is  not  limited  by  the  thickness  of  the  polymer  layer  since  it  is 
measuring  lateral  diffusion  as  opposed  to  through-thickness 
diffusion.  Additionally,  the  time-scale  of  the  test  is  determined  by 
the  contact  radius,  which  is  experimentally  controlled  by  the 
indentation  depth.  Therefore,  coupling  with  other  relaxation 
processes  can  be  avoided  because  the  time-scale  for  poroelastic 
relaxation  can  be  adjusted.  Since  PRI  testing  couples  contact 
mechanics  to  fluid  transport,  the  indentation  profile  will  deter¬ 
mine  the  migration  behavior  of  the  solvent.  Thus,  the  key 
assumption  for  extracting  the  materials  properties  from  PRI  is 
that  the  mechanical  and  transport  properties  of  the  polymer  layer 
of  interest  are  considered  to  be  isotropic. 

Conclusions 

We  have  demonstrated  PRI  as  an  alternative  and  simple  testing 
approach  for  characterizing  the  transport  properties  of  thin  PEG 
hydrogel  layers.  By  utilizing  a  linear  poroelastic  theory  that 
predicts  the  relaxation  response  over  a  range  of  geometric 
confinement  from  unconfined  to  highly  confined,  the  results  from 
PRI  testing  enable  determination  of  three  materials  properties  of 
the  hydrogel  including  shear  modulus,  diffusion  coefficient,  and 
average  pore  dimensions. 

In  the  highly  geometrically  confined  scenario,  the  PRI 
approach  is  quantifying  the  lateral  diffusion  process.  In  theory, 
PRI  is  not  limited  by  the  thickness  of  the  polymer  layer.  The 
important  aspects  to  consider  are:  (1)  the  materials  properties  ( D , 
G)  of  the  polymer  are  assumed  to  be  isotropic,  and  (2)  the  contact 
radius  is  sufficient  large  to  avoid  overlapping  with  other  relax¬ 
ation  time-scales,  which  can  be  accomplished  via  increases  in  R 
or  5.  While  this  work  has  demonstrated  that  PRI  can  measure  the 
transport  for  “thin”  or  geometrically  confined  hydrogel  layers, 
more  work  is  needed  to  show  that  the  technique  can  measure 
thinner  polymer  layers  such  as  membranes  used  in  water  purifi¬ 
cation.  Additionally,  we  are  interested  in  imaging  the  deformed 
region  of  the  hydrogel  during  PRI  testing  to  validate  the 
poroelastic  process  as  a  function  of  geometric  confinement. 

Another  point  we  would  like  to  make  is  that  many  polymers 
gels  behave  as  hyperelastic  materials  at  large  strains.  Thus, 
a  systematic  study  of  the  effect  of  nonlinear  strains  on  their 
poroelastic  relaxation  process  is  warranted.  Work  is  underway  to 
study  the  effects  of  nonlinear  strains  on  the  poroelastic  relaxation 
of  physical  gels  such  as  block  copolymer  photonic  gels  where 
their  mechanochromic  properties  are  strongly  tied  to 
poromechanics . 39 

Finally,  based  on  eqn  (1),  the  poroelastic  relaxation  time  is 
linked  with  the  diffusion  coefficient  of  small  molecules.  There¬ 
fore,  PRI  should  be  capable  of  measuring  the  diffusion  of 
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multiple  small  molecule  species  within  the  polymer  layer, 
provided  that  the  diffusion  coefficients  are  unique  and  the 
diffusion  processes  are  uncoupled.  This  measurement  can 
potentially  be  applicable  for  quantifying  the  diffusion  of  phar¬ 
maceuticals  within  hydrogel-based  drug  particles.  Work  is 
currently  underway  to  characterize  the  diffusion  of  water  and 
linear  polymer  chains  within  hydrogel  layers  using  PRI. 
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